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Abstract 

Following recent studies of non-relativistic reductions of the single-nucleon elec- 
tromagnetic current operator, here we extend the treatment to include meson ex- 
change current operators. We focus on one-particle emission electronuclear reac- 
tions. In contrast to the traditional scheme where approximations are made for the 
transferred momentum, transferred energy and momenta of the initial-state struck 
nucleons, we treat the problem exactly for the transferred energy and momentum, 
thus obtaining new current operators which retain important aspects of relativity 
not taken into account in the traditional non-relativistic reductions. We calcu- 
late the matrix elements of our current operators between the Fermi sphere and 
a particle-hole state for several choices of kinematics. We present a comparison 
between our results using approximate current operators and those obtained using 
the fully-relativistic operators, as well as with results obtained using the traditional 
non-relativistic current operators. 

PACS numbers: 25.20.Lj, 25.30.Fj, 25.10+s 
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1. Introduction 



In recent work |T], 0, |3|] an improved version of the single-nucleon electromagnetic current 
has been studied. There the so-called "on-shell form" of the current was derived as a non- 
relativistic expansion in terms of the dimensionless parameter rj = p/m^, where p is the 
three-momentum of the struck nucleon (the one in the initial nuclear state to which the 
virtual photon in electron scattering reactions is attached) and is the nucleon mass. 
Generally (that is, for nucleons in typical initial-state nuclear wave functions) p lies below 
a few hundred MeV/c and thus rj is characteristically of order 1/4. Barring some extreme 
choice of kinematics such as the selection of extremely large missing momenta in (e, e'N) 
reactions — conditions for which presently no approach can be guaranteed to work - 
an expansion in powers of rj is well motivated. Similar arguments do not however apply 
equally for other dimensionless scales in the problem. Indeed, a specific goal of this past 
work has been to obtain current operators which are not expanded in either k = q/2rriN 
or A = u/2m^, where q is the three-momentum and to the energy transferred in the 
scattering process, since one wishes the formalism to be applicable at GeV energies where 
these dimensionless variables are clearly not small. 

Traditionally, many studies have indeed been undertaken assuming that k < 1 and 
A < 1 aimed of course at treatments where non-relativistic wave functions are em- 
ployed F° r high-energy conditions the current operators so obtained are bound to 
fail, whereas our past work on the single-nucleon current provides a way to incorporate 
classes of relativistic corrections into improved, effective operators for use with the same 
non-relativistic wave functions. 

Not only the single-nucleon (one-body) current, but also the two-body meson ex- 
change currents (MEC) have frequently been evaluated using similar traditional non- 
relativistic expansions ||-[21j in which k and A are both treated as being small, together 
with the assumptions that all nucleon three-momenta in the problem are small compared 
with m N . In other work f22||-p5[l, relativistic currents have been used directly in cases 



where the nuclear modeling permitted. 

Our goal in the present work is to extend our previous approach for the single-nucleon 
current operators now to include a treatment of pion-exchange MEC. We make expansions 
only in rji = pi/m^, where {p{\ are the initial-state nucleon three-momenta, whereas 
we treat the dependences in the on-shell form exactly for k, A and any high-energy 
nucleon momenta, specifically for any nucleons in the final state not restricted to lie 
within the Fermi sea. Such new MEC operators may straightforwardly be employed in 
place of previous non-relativistic expansions using the same non-relativistic initial and 
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final nuclear wave functions employed in the past, since our effective operators incorporate 
specific classes of relativistic effects (see the discussions of the single-nucleon current 
referred to above). 

In the present work, as a first step, we focus on the general form of the MEC matrix 
elements for pionic diagrams (the so-called seagull and pion-in-flight contributions) and 
plan to extend our treatment to other diagrams and other meson exchanges in future work. 
We do not present any results for electromagnetic response functions, postponing such 
discussions until the corresponding correlation effects have been brought under control 
(also work in progress). Finally, in the present work we focus on specific classes of matrix 
elements, namely those with one high-energy nucleon in the final state, i.e. one-particle- 
one-hole (lp-lh) matrix elements; in subsequent work we shall extend the scope to include 
2p-2h configurations. 

The organization of the paper is as follows: After reviewing the treatment of the 
single-nucleon current in Sect. 2.1, in Sect. 2.2 we discuss the new approach for the elec- 
tromagnetic meson-exchange currents that treats the problem exactly for the transferred 
energy and transferred momentum. We check the quality of our expansions in powers of 
the bound nucleon momenta divided by by calculating the matrix elements of the 
current operators between the Fermi sphere and a particle-hole state. We compare with 
the matrix elements obtained using the full current operators as well as with the results 
for the traditional non-relativistic expansions. These results are presented in Sect. 3 
(with reference to specific details that are covered in an appendix), together with a brief 
discussion of the high-g limits reached by the currents. Finally in Sect. 4 we summarize 
our main conclusions. 



2. Current Operators 

2.1. The Electromagnetic Current Operator 

We start our discussion with the single-nucleon on-shell electromagnetic current oper- 
ator and its non-relativistic reduction. Although this case has been already treated in 
detail in Refs. [I], 0, |3|, here we provide a simpler derivation that will be also applied to 
the case of the MEC operators for one-particle emission reactions. The single-nucleon 
electromagnetic current reads 



J"(PV; Ps) = u(p', s' 



Fi{Q 2 )r + ~ F 2 {Q 2 )a^Q v 



2m 



N 



u[p, s) 



(1) 
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where P M = (E, p) is the four-momentum of the incident nucleon, P' M = (E', p') the 
four-momentum of the outgoing nucleon and = P' M — P^ = (u, q) the transferred 
four-momentum. The spin projections for incoming and outgoing nucleons are labeled s 
and s', respectively. We follow the conventions of Bjorken and Drell |26| for the w-spinors. 



For convenience in the discussions that follow of the scales in the problem we introduce 
the dimensionless variables: r\ = p/m N , e = E/m N = yl + rf, A = i^^-, k = and 
t = — 7^2- = k? — A 2 . For the outgoing nucleon, rj' and e' are defined correspondingly. 



iY 



For any general operator whose 7-matrix form is given by 

r-^ 11 i"\ (2) 



one has u(p', s')Tu(p, s) = xl'^Xs, with the current operator T given by 

— 1 r, T, T / <t ■ ri a ■ ri' <r ■ ri 1 cr ■ ri\ 

T = -J 1+5 fl + £' ) Tn + Tia -T21 -T22 ' • 3 

2 VV A ; V 11 12 1 + e 1 + e' 1 + e' 1 + eJ y 1 

This general result will be used throughout this work in discussing the non-relativistic 
reductions of the various current operators. 

An important point in our approach is that we expand only in powers of the bound 
nucleon momentum 77, not in the transferred momentum k or the transferred energy 
A. This is a very reasonable approximation as the momentum of the initial nucleon is 
relatively low in most cases, since the typical values of 77 lie below r\p = kp/m^, where kp is 
the Fermi momentum (77^ is typically about 1/4). However, for those cases corresponding 
to short-range properties of the nuclear wave functions it will be necessary to be very 
careful with the approximations made. Indeed, for large values of 77 a fully-relativistic 
approach will likely prove necessary. Expanding up to first order in powers of 77 we get 
e ~ 1 and e' ~ 1 + 2A. Thus, the non-relativistic reductions of the time and space 
components of the single-nucleon electromagnetic current operator can be evaluated in a 
rather simple form. 

Let us consider first the case of the time component. We have 

J°(P's'; Ps) = H(p', s')J°u(p, s) = xl>J*Xs , (4) 

with the current operator J° = P17 + iF2<J 0u Q u /2rriN- Using the general result given by 
Eq. (^j and expanding up to first order in 77, it is straightforward to get the relation 

J* ~ -^G E + (g m -£f)( KX f,).tr, (5) 
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where we have introduced the Sachs form factors Ge = F\ — tF 2 and Gm = F\ + F 2 , and 
have used the relations 



A ~ r + n ■ rj (6) 
k 2 ~ r(l + r + 2k • 77) . (7) 

The expression (|5]) coincides with the leading-order expressions already obtained in pre- 
vious work [0, 0; in those studies a different approach was taken which, while more 
cumbersome, does yield terms of higher order than the ones considered in the present 
work. It is important to remark again that no expansions have been made in terms of 
the transferred energy and transferred momentum; indeed, K, A and r may be arbitrarily 
large in our approach. 

Let us consider now the case of space components. Thus, we have J(P's'; Ps) = 
u(p', s')Ju(p, s) = xl'JXs- Introducing the matrix form of the vector component for the 
single-nucleon electromagnetic current operator in the general relation ([3]), one can finally 
write 

iG M (cr x k) + ( G E + -G M ) rj + G e k 



vT+7 1 ' \ " 2 

Gm 1 \ ^Ge / \ 

(K, ■ T))H, -{a X K,)K, ■ 77 



2(l + r) v " 2(1 + r) 
- ir(G M - G E /2){* x 77) + Z( ^~^ (K x 77) CT ■ k\ , (8) 

where we have used the relations given by Eqs. (||,0). 

In order to compare with the previous work 0, we write the expression for the trans- 
verse component of the current, i.e., J" 1 " = J — ^rK>. After some algebra we get the final 
result 




(9) 

It is straightforward to prove that this expression coincides with the result given by 
Eq. (25) in Ref. pj for an expansion in powers of 77 up to first order. 

Therefore, as can be seen from Eqs. (§,0), at linear order in 7] we retain the spin-orbit 
part of the charge and one of the relativistic corrections to the transverse current, the 
first-order convective spin-orbit term. It is also important to remark here that the current 
operators given by Eqs. (|5|||) satisfy the property of current conservation AJ = n ■ J . 
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Finally, it is also interesting to quote the results obtained in the traditional non-relativistic 
reduction 0,0-0, where it is assumed that k « 1 and A << 1: 



nonrel Gjj] 
~J nonrel = -iG M [n X a] + G 1 







J]- 









(10) 



Note that this traditional non-relativistic reduction contains both terms of zeroth and 
first order in t], i.e., the convection current, and is therefore not actually of lowest order 
in r). 

The present expansion for the electromagnetic current operator of the nucleon was 
first checked in Ref. |l|], where the inclusive longitudinal and transverse responses of a 
non-relativistic Fermi gas were found to agree with the exact relativistic result within a 
few percent if one uses relativistic kinematics when computing the energy of the ejected 
nucleon. Recently the same expansion has been tested with great success by comparing 
with the relativistic exclusive polarized responses for the 2 H(e, e'p) reaction at high mo- 
mentum transfers ||. This relativized current has also been applied to the calculation of 
inclusive and exclusive responses that arise in the scattering of polarized electrons from 



polarized nuclei 0. 27 



We see that the expansion of the current to first order in the variable r] = p/m^ yields 
quite simple expressions; moreover the various surviving pieces of the relativized current 
(i.e., charge and spin-orbit in the longitudinal and magnetization and convection in the 
transverse) differ from the traditional non-relativistic expressions only by multiplicative 
(g, ^-dependent factors such as k/^/t or 1/vT+T, and therefore are easy to implement 
in already existing non-relativistic models. In the next section we perform a similar 
expansion for the MEC and later return to check our results through direct comparisons 
with the exact relativistic matrix elements. 



2.2. Meson— Exchange Currents 

Once the procedures for expanding the single-nucleon electromagnetic current are fixed, it 
is clear how to proceed to obtain relativistic expansions for the meson-exchange currents. 
In the present work we begin by focusing on pion exchange MEC effects, leaving the 
treatment of other mesons to future work. Following the ideas and methods developed 
in the previous section, our main aim here is to get new non-relativistic reductions for 
MEC treating the problem of the transferred energy and transferred momentum as above, 
namely in an un-expanded form while expanding only in the initial nucleon momenta. In 
this way the expressions obtained will retain important aspects of relativity not included 



6 



in the traditional non-relativistic MEC used throughout the literature. 

Let us consider the meson exchange current operator Jff . For defmiteness we focus 
on one-particle emission reactions where the matrix element of Jff EC taken between the 
Fermi sphere and a particle—hole state, namely 

(ph-vr c \ F )=T, e [( P hvr c \ hh ') - (ph'\j? Ec \h f h)] , (ii) 

s',f h'<k F 

kp being the Fermi momentum, is the relevant one. Since, as seen below when the pion- 
exchange operators are given, the currents being considered all have isospin dependences 
of the form r a (1)t&(2), the first term (direct term) vanishes, 

J2(t P t'\r a (l)r b (2)\t h t') = (t p \r a \t h ) J2(t'\r b \t') = . (12) 

v v 

Therefore, the only remaining term is the exchange term, and we can simply write 

{ph ^\jMEC lF) = ^ <ph>\J?BC\h'h) . (13) 

s',t' h'<k F 

In what follows we will be interested in the evaluation of the particle-hole matrix elements 
{ph'\J^ IEC \h'h) and their new non-relativistic expressions. As in our treatment of the 
single-nucleon current, it is convenient to express the results in terms of spin matrix 
elements of particular operators : 

{ P h'\j s /\h'h) = x\xl' 2 Jl' P Mx sl x S2 , (14) 

where S(P) denotes the seagull (pion-in-flight) contributions to the MEC, as shown in 
Fig. 1. 

2.2.1. Seagull Current Operator 

The relativistic seagull current operator is given by 

J »(Q) = y2~2 iezabU (P' 'n s 'i) T al5 fiMpi, si) R2 j m2 u(p' 2 , s' 2 )r b ^ lx u(p 2 , s 2 ) 

+ (1 »2) , * (15) 

where V is the volume enclosing the system and the different nucleon kinematic variables 
are: P l = (E 1 ,p 1 ), P{ = P 2 = (E 2 ,p 2 ) and P' 2 = (E 2 ,p' 2 ) (see Fig. 1). The 

four-momentum of the exchanged pion is K x = (E kl , k x ) (with K 2 likewise) and its mass 
is m n . The terms / and represent the pion-nucleon coupling and pseudovector form 
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factor, respectively. Note that the following kinematic relations are satisfied for the two 
diagrams involved 



rfi 
r?' 2 



Ci 

C-2 



r\ 2 + 2k 
Vi 

v'i - Vi = V2 - m + 2k 

V 2 - V2 = Vl - f)2 , 



(16) 
(17) 
(18) 
(19) 



where we follow the general notation introduced in Section 2.1 and have also defined 
£ 12 = ki i2 /m,jv with k 12 the three-momenta of the exchanged pions. 

The particle-hole matrix element of the seagull current is then given by 

P 



(ph'Mh'h) 



V 2 m; 



:i£zab(t p \T a \th>)(th'\n\th) 



U (~P'2, 4)75 -^2«(P2, S 2 ) 



v 



F? 



Ill; 



-u(p'2i S 2 h^MP2' S 2) 



K 2 — rri; 



(20) 



where now p[ = p is the momentum of the ejected particle above the Fermi sea, which 
can be large for large momentum transfer, while P2 = h is the momentum of the bound 
nucleon before the interaction (related to the missing momentum in (e,e'p) reactions), 
which can be considered small compared to the nucleon mass. Finally, p 2 = Pi = h' is 
the intermediate momentum of the bound nucleon interacting with the ejected nucleon 
by pion exchange. Therefore p 2 and pi are small compared with the nucleon mass. In 
dimensionless terms, we can safely expand in rf 1 , r} 2 and hence r]' 2 , whereas we cannot in 
general expand in r][. Moreover, £ 2 is small, whereas £i can be large. By analogy with 
the single-nucleon electromagnetic current, let us now proceed to the evaluation of the 
time and space components separately. 

Time Component 

Using the general result given by Eq. (|) and the matrices 75 fii and 7570 we can write 

f 



x 



X 




AV 2 ml 



ie Z ab{t P \T a \t h ,){t h ,\T b \t h ) ^/(l + £!)(! + e 2 )(l + ei)(l + e' 2 ) 



m N a- ■ Ci - ^k 1 — — - ^k 1 — -7 + m N ——{a- ■ Ci)t— — 

1 + E\ 1 + £7 1 -r £j 1 + Si 



a- ■ rfj 



a ■ 7)1 



Xsi 



& ■ *)2 , g • V2 

1 + E 2 l + e' 2 



X 



S2 



8 



X 



X. 



1 + ei l + e[ 



Xsx 



v 



mi 



&2 



12 



l+e 2 



E k 2 — 7 + m N 



1 + 4 



1+4 



& ■ Q 2 



V 2 



1 + e 2 



Xs 



(21) 



where E^E^) is, as mentioned previously, the energy of the pion. 

We note that the resulting expansion for the MEC should be used together with 
the single-nucleon current, which has been developed to first order in rj. Therefore, 
in order to be consistent, we want to perform the expansion of the MEC also to first 
order in the corresponding small quantities {ry 1 , r] 2 , rj 2 , £ 2 }, whereas {rj^, k, £ x } are treated 
exactly. Using the kinematic relations given in Eqs. (p^6|-p!9|) , the following non-relativistic 
kinematic reductions are involved (up to first order in the small quantities) 



e 2 ~ e 2 ~ 



e[ ~ 1 + 2A 
E kl ~ 2m N A 
E k2 ^ 0, 



(22) 
(23) 
(24) 
(25) 



where we follow the notation introduced in Section 2.1. Using these non-relativistic 
expansions and the kinematic relations given by Eqs. (|iT)|-|iTS|) , it is straightforward to 
obtain 

F f o-i ■ (Ci - A*7i) cr 2 ■ (t/ 2 + rf 2 ) 



4(1,2) 



2VTTX 1 K\ - m\ 

CTf [(1 + A)T|! + Tfi] CT 2 -C 2 \ 



mi 



having defined 



f 2 m N 
V 2 mt 



{t P \r a \t h '){t h '\r b \t h )F 1 . 



(26) 



(27) 



We can still simplify this expression by using the relation A ~ r + k, ■ rj 2 (valid up to first 
order in powers of rj 2 ): the result is 



4(1,2) 



F \ °"i ■ (Ci - ™7i) o-2 ■ {l 2 + V 2 ) 



2VT 



mi 



CTi ■ [(1 + T)T] X + T7i] CT 2 -C 2 



mi 



(28) 



Examining this result, we see that we have retained some terms that should actually be 
neglected because they give 2 nd -order contributions: this is the case for the factor t(<ti-?7 1 ) 
which is multiplied by a 2 ■ (r) 2 + 772) and a 2 ■ £ 2 . When these contributions to the current 
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are omitted, then the comparison between our result and the traditional non-relativistic 
expression |I5|| , [ |TTf becomes more straightforward. Moreover note that if one neglects 
the term t(ct 1 ■ T7-J, then one recovers for the time component of the seagull current an 



expression that is similar to the traditional non-relativistic reduction ||T5| except for the 



common factor 1 / y/1 + t, which accordingly incorporates important aspects of relativity 
not considered in the traditional non-relativistic reduction. This result is similar to the 
discussion given in Ref. Jl[j for the case of the single-nucleon electromagnetic current. 
Finally, note that, strictly speaking, in order to be consistent with the non-relativistic 
expansion, the contribution of (j\ ■ rf l should be also neglected. Therefore, our final 
expression for the non-relativistic reduction of the seagull current can be written in the 
form 



4(1,2) 



(a-x ■ Ci) <r 2 • (t7 2 



Oi • y/i) O2 • C2) 



Kl - ml 



2^1+^ { K( - ml xv 2 

In order to obtain a truly first-order expansion of the current it is convenient to re- 
express the momenta involved in Eq. (|28|) in terms of the momentum transfer k, which 
can in principle be large, and the nucleon momenta rj 1 , r\ 2l which lie below the Fermi 
surface and are kept as the parameters of the expansion. By inserting Eqs. (|16|-[T9|) into 
Eq. (|28|) and keeping only terms linear in rf^ rj 2 , one gets (note that we keep the full pion 
propagators) 



7*0,2) = 



<r 2 • (77i + il 2 ) cr 2 ■ (r] 1 - r] 2 ) 



mi 



mi 



(30) 



Finally, it is also interesting to examine the limit 77^ — 0, since this will provide some 
understanding of how the MEC effects are expected to evolve in going from light (rjp very 
small) to heavy nuclei (t)f = 0.29). Obviously, in this case the following relations are 
satisfied 



Vi = r} 2 = r) 2 
Ci = m = In 



and therefore the seagull current simply reduces to 







4(1,2) 



. 



(31) 
(32) 



(33) 



This is a consequence of the fact that the time component of the seagull current is of first 
order in the small variables involved or, in other words, it is of 0(t)f)- 



Space Components 
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The particle-hole matrix element is given by 

P 



(ph'\J b \h'h) 



V 2 m: 



-itzab (tp\r a \t h ') {t h > \r h \t h ) 



x S«(p'i, s'i)T5 -Kl«(Pl,Si) 



w(p'2> 4)757M(P2, s 2 ) 



- "(p'a, 4)75 fe2u(p 2 , s 2 )^p r 1 — ^(P'li s 'l)757^(Pl, «l) 



' Kl - ml 



(34) 



Using again the general relation (12) for the matrix forms for 757 and 75 fii } 2, one can 
write 



(ph'\J s \h'h) 
P 



4V 2 mi 



itzab (t P \T a \t h ,)(t h ,\n\t h ) v /(l + ei)(l+e 2 )(l+ e[)(l + 4) 



cr ■ 77 x cr • 



Si 



X 



if 



- X 



K\ - ml 
t 



X 



cr H -<x 



1 + ei 1 + £1 



An — ml 



CT-r7 2 o-77 2 



miv<T ■ C 2 - E k2 — — + — — + m N ——cr ■ £ 2 - 

V 1 + £ 2 1 + eW 1 + ^2 1 



We can now make use of the non-relativistic reductions as explained in the previous 
section in connection with the time component. Then, considering only terms up to first 
order in the small quantities in which we are expanding, after some algebra one finally 
gets 



J (1,2) ~T- 



1 o-a • (Ci - At/0 <x 2 



vTTa 



°"l ( <T 2 • C2 

Ko — mi 



(36) 



^v/TTA Kl-ml ,v 2 
This expression can be further simplified by using the relation A ~ t+k-77 2 and expanding 
the terms 1/yl-j-A and y/1 + A in powers of 772 (up to first order). The final result is 



J (1,2) ~T 



1 (Ti ■ 



1 " 2775) ) Cl - T7) 



2(l+r) 



<T2 



Vl+7 



K\ - ml 



(gj • C2 

K| - ml 



(37) 



It is important to note that neglecting the term (k • 77 2 )/[2(l + r)] compared to 1, and 
r(cr 1 ■ r/ 1 ) compared to <j\ ■ £ x (good approximations — see the next section) one simply 
recovers the traditional non-relativistic expression |T3| except for the factors 1/y/l + r 
and \/l + t that multiply the contributions given by the two diagrams involved. As in the 
case of the time component, this result indicates that important relativistic effects can 
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be simply accounted for by these multiplicative terms. In next section we shall present 
results for a wide choice of kinematics showing the validity of the obtained expressions. 
By inserting Eqs. ([T6|-|l9|) into Eq. (^) and expanding up to first order in t] 1 , rj 2 , one gets 

2(Ti ■ K 



J S (h2) 



01 



1 - 



K ■ T) 2 



+ 



<Ti • rj 2 



2(1 + r) v/TT7 



vT 



rax ■ rj 1 



02 



K? — m 2 

1 7T 



Finally, in the limit rjp — > 0, one obtains 



j s hx 



IT 



,<Ti ■ K)<7 2 



(38) 



(39) 



»?f^o v /TTr Q 2 -ml ' 

which shows that the space components of the seagull current are of 0(1) and contribute 
even for nucleons at rest, as do the charge and magnetization pieces of the one-body 
current. 



2.2.2. Pion— in— flight Current Operator 

The relativistic pion-in-flight current operator reads 



JP ^ Q) = V^ ezah \K\ - ml)(K$ - ml) 

where the kinematic variables are defined in Fig. 1 and where the kinematic relationships 
given m Eqs. (HMD for the seagull diagram are again satisfied. In order to preserve 
gauge invar iance, we choose the electromagnetic pion form factor to be F w = FY (see 
Ref. 0). 

The exchange particle-hole matrix element is given by 



{K x - K 2 )^ 



-ufa'-^TaryB ^lu(Pl)u(P2) T bl5 ^2u(p 2 ) , (40) 



P F v (K 
(ph'\J^\h'h) = ~ j^-^ ie zab (t p \ r a \ t h > )(t h >\r b \t h ) 



K. 



2)n 



V 2 ml 



(K 2 - ml){K\ - ml) 



X w(Pl)75 ^im(Pi)m(P 2 )75 ^2«(P2) • 

Again using Eq. (fj) and the matrix forms of 75 fei (75 ^2) we can write 



(41) 



p F v 

(pti\J^\tih) = A T r9 _ 9 ie zab {t p \r a \ t h , ){t h >\n\t h ) 



AV 2 ml 



x 



{K, - K 



2Jm 



x X 



(1 + £l )(l + e[)(l + s 2 )(l + ^ {K 2_ m 2 ){K 2_ ml) 

' rr t F ( °" Tyi 4- ±m ° ' ^ 

m N er ■ Q x — h, kx \ — h — — - J + mjv 



x Xs 



m N a- ■ C 2 ~ E, 



k 2 



° n 2 , °" »72 



rT7[ (<T ' Cl) rT^ 



■12 I ^ ' 'H2 I /• , " 

TT^ J m7V TT^ ' 1 



• *?2 



1 + 6 2 1 + eS 



£2 



Xsi 
Xs2 



(42) 
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Following the arguments discussed for the seagull current we expand up to first order 
in powers of the variables {rf^ r] 2 , rj 2 , C 2 }> whereas {Vij^jCi} are treated exactly. The 
non-relativistic reductions for the kinematic variables are given by Eqs. p^-p5|) and after 
some algebra one obtains 

J ' (1 ' !) = "7BMM' rtl ' rC! ' (43) 



Using again the first-order relation A ~ r + k, ■ rj 2 and expanding I/a/1 + A in powers of 
r]2, the p ion-in-flight matrix element may finally be cast into the form 

-jP n x F m N (K x - K 2 ) li . . 



This expression is similar to the traditional non-relativistic current |T5] except for the 
common factor 1 / y/1 + r, which should again include important aspects of relativity not 
taken into account in the traditional non-relativistic reduction. 

Once more we can express this matrix element in terms of k, t] 1 and r] 2 and keep only 
linear terms in the small momenta, obtaining 

—p. 4J-m% r 

J ° (1 ' 2) = " Wl - mi) {Kl - ml) VTTT ' 71 • ^ • (?h " ^ (45) 
j p (l 2) - <n ■ ^2 ■ fa - r? 2 ) 

Note that the space component of the pionic current is, in leading order, purely longitu- 
dinal; its transverse components are in fact of second order in rjp. 

Finally the limit tjf — > reduces to (J,f) = 0, and we see that all the components of 
the pionic current are of 0{t]f) m the expansion. 



3. Results 

We present here a discussion of the numerical results obtained for the MEC matrix ele- 
ments. In the following we compare the fully-relativistic results with those obtained using 
two different expansions for the MEC: i) the traditional non-relativistic approach (TNR), 
where k, A and all nucleon three-momenta (rj) are treated as being small and ii) our new 
non-relativistic (also referred to as "relativized") approach (NR) where we expand only 
in powers of bound nucleon three-momenta, not in k or A. In order to check the validity 
of our expansions we compute the transition matrix element of the current between the 
Fermi sea and a p-h excitation, i.e. {phr 1 \jff EC \F). Furthermore, to assess the quality 



13 



of the expansions performed for the MEC in the last section, we show here comparisons 
of matrix elements for the relativistic MEC taken between Dirac spinors with matrix el- 
ements of the expanded MEC taken between Pauli spinors. In what follows we use the 
notation of the lower diagrams shown in Fig. 1. The connection with the general termi- 
nology introduced in the previous section is given by the following relations: P = P[ is 
the four-momentum of the ejected nucleon and H = P 2 , H' = P' 2 = Pi are the initial and 
intermediate four-momenta of the bound nucleons. The corresponding three-momenta 
are denoted by: p = |p| = \p[\, h — |h| = | p>2 1 and h' = |h'| = | p>2 1 = |pi|- 



3.1. Seagull current 



In order to simplify our analysis we first extract from the currents the factors which are 
common to both relativistic and non-relativistic currents, namely, coupling constants, 
form factors and isospin matrix elements. Accordingly, for the seagull current J~, we 
define a dimensionless function K^(q, u,h) as follows: 

P 



T,{pti\J s Jh'h) 



Vml 



i([TiXT 2 ] z )F?k? F K?(q,u,h) 



(47) 



where ([ti x T2\ z ) stands for the corresponding matrix element and implies a summation 
over isospin. Using the definition of the seagull current, the expression for the function 



d 3 h' 
(2k) 



-y 



u(ti)^($'- #)u{h)_ 



u(p)j 5 (P- fl')u{ti) 
(P - H'Y - ml 



(H' - H) 2 - ml 



U{v)lhly,u{ti) 



(4£ 



where the sum runs over the third spin component of the spinor u(h') = u(h!, s^/) and 
the integral over h' is performed below the Fermi sea |h'| < hp. This expression for is 
the exact relativistic result. Note that we have divided by the squared Fermi momentum 
kp in order to obtain a dimensionless function; a factor k F is correspondingly included in 



Eq. (i7[). In the previous section we have performed an expansion in powers of the small 
quantities h/m^ and h'/m^- Therefore we define for the components of the function 
the following "non-relativistic" approximations : 

d 3 ti 



2m N k 2 F Vl + r 
XL* • (p - h' - 



E 



(2k) 



rh') Xsh , 



(P-H 



A2 



mi 



(h' + h) Xsh 
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fh' 



K s > Nm (q,u,h) = 



(H' - H) 2 
1 1 



• (P + h' + rh') Xshl 



mi 



(49) 



E 



• d 3 ti 



(» 



(2tt) 3 

K-h 



2mjv(l+r) 



fp - h') - rh' 



-a 



(P - P') 2 - "4 

(h' ~ h) Xsh f 

2 -^Sp* 7 -^^' 



mi 



(50) 



where "S,NR1" is meant to denote non-relativistic approximation number 1 for the seagull 
contributions. These expressions should be compared with the traditional non-relativistic 
seagull current that can be obtained by taking the limit k — > and r — > 0, namely 

d 3 h' \xl p <r ■ (P ~ h!) X s h , f 



jy-S,TNR/ i \ 



2m ^k F 



E 



fh' 



(27T) 



(if' - if) 2 



(P - P') 2 

xLo" • (p + h ')x Sh , 



(51) 



K 



S.TNR 



1 



7.2 E 



d 3 /i' rxi p ^ ■ ( P - hQx^ t 



(2k) 



(P - H') 2 - 

xl h °- (h'-h)xs h f 
(P'-P) 2 -m 2 Xs p aX ^' 



mi 



Xs h , (T Xs h 



(52) 



where "S,TNR" denotes the traditional non-relativistic approximation for the seagull 
contributions. In addition one would like to find approximations to the currents where 
the relativistic effects are accounted for as corrective factors consisting of simple functions 
of (q, uj) such as the combination of k, t and y/1 + r previously found for the case of the 
single-nucleon current (see Sect. 2.1). One could thus easily implement the relativistic 
corrections in existing models of traditional non-relativistic MEC. For these reasons we 
define a second approximation for the seagull current in which we neglect the factors n ■ 
h/[2mjv(l + r)] and rh' in the function K ,NRl , thereby yielding a second non-relativistic 
approximation for the seagull current, K^' NR2 , which differs from the traditional current 
only by the factor y/1 + r: 

d 3 /V I xjo- ■ (P ~ h') Xa 
2m N k 2 F y/T^T tT J (2vr) 3 [ (P - H' f 



K, 



S,NR2 




mi 



V + 
—Xl.v 



(h' + h) Xsh 



Xl h ,<r ■ fa' - h)xs h f 



(P'-P) 2 -m 2 ASp 



xLo- ■ (p + h') Xt 



h> 



1 



y/l + r 



S.TNR, 



q,u,h) 



(53) 
(54) 
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K S ' NR2 (q,u,h) 



1 1 



+ 




2 _ m 2 Xs^^Xsh 



'7T 



(55) 



where "S,NR2" denotes non-relativistic approximation number 2 for the seagull contribu- 
tions. In what follows we check the validity of the various approximations to the full cur- 
rent introduced above by performing numerical calculations of the functions K^(q,u,h) 
for several choices of the kinematical variables. 

First we notice that, for fixed momentum and energy transfer, (q,u), there are re- 
strictions on the values of the momentum of the hole h, since our nucleons are on-shell 
and the momentum of the ejected particle P M = H* 1 + must satisfy P^P^ = m 2 N . In 
particular, the restriction 



fixing the angle between the hole momentum h and the momentum transfer q is seen to 



is the magnitude of the hole momentum and 4>h is the azimuthal angle of h in a coordinate 
system with the z-axis in the direction of q. The angle between h and q is then given by 



As this must lie between —1 and 1, one obtains a restriction on the values of uj for which 
a contribution to the on-shell matrix element exists. 

For fixed values of (q,u), there is another restriction generated by the condition that 
the particle momentum p must lie above the Fermi sea; in this work, however, we are 
only interested in the high-momentum region where relativistic corrections are expected 
to be important and there such Pauli-blocking effects can be ignored. 

In Figs. 2-4 we show the dominant (real or imaginary) parts of the four vector com- 
ponents of the seagull function K^. The components not shown in the figures are found 
to be negligible in our calculations, as a result of cancellations occurring among differ- 
ent pieces: in order to understand the reasons for these cancellations we have explored 
the symmetries of the integrals involved in the various components. A summary of that 
study is given in the Appendix for the particular case of the pion-in-flight current; a 
similar procedure can be followed for the seagull current. 

We choose the typical value fop = 250 MeV/c for the Fermi momentum and in all of 
the figures the hole kinematics correspond to h = 175 MeV/c and 4> h = 0° . Because of the 



2h • q = uj 2 - q 2 + 2E h u , 



(56) 



hold. Therefore the functions depend only on the variables (q, uj, h, (f) h ), where h — |h 



cos6^ = 



to 2 - q 2 + 2E h uj 
2hq 



(57) 
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above mentioned symmetries in the currents it is possible to relate the matrix elements 
corresponding to other choices of the angle <ph to the ones calculated here — for example 
the symmetry between the <ph = 90° and fa = 0° cases is discussed for the pion-in-flight 
current in the Appendix. We have checked that the same connection obtains for the values 
<ph = 90° and <ph = 180°, with the exception that in this case the roles of some of the 
components are switched. However, since no new information concerning the validity of 
the expansion is obtained from values of fa different from zero, we show results only for 
the latter case. 

The function is displayed for three values of the momentum transfer, namely 
q = 500, 1000 and 2000 MeV/c, as a function of u. For each q, the allowed values of uo are 
restricted to the intervals displayed in the figures. Note that the cj-values in the figures 
have been chosen to lie in a region around the approximate quasielastic peak position, 



oo = y q 2 + m 2 N — m^, which for the selected momenta occurs for lo ~ 125,433 and 
1271 MeV, respectively (approximately at the center of the u;-region displayed in each 
situation). 

In each of the panels in Figs. 2-4 we plot four curves, corresponding to the fully- 
relativistic (solid lines) and to the non-relativistic approximation K^' NRl (dashed 
lines); moreover the results obtained with the traditional non-relativistic current K~> TNR 
(dot-dashed lines) and with our simplified, non-relativistic approximation K^' NR2 (dotted 
lines) are also shown. The functions are spin-matrices, i.e., they depend on the spin 
projections s p , Sh of the particle and of the hole, respectively; accordingly we write 



For the sake of brevity, in the figures we show results only for the spin components K\\ 
and K\2, corresponding to (s p ,Sh) = (1/2,1/2) and (1/2,-1/2) respectively. Similar 
results are found for the remaining components of the currents. As for the dependence 
upon the angles fa, relationships between the spin components displayed in the figures 
and the remaining ones can be established and the same comments made with respect to 
the dependence on fa are valid here as well. 

Looking at Figs. 2-4, we first note that the imaginary parts of Kq , Kf and K% and the 
real part of iff are not shown. Indeed, as mentioned before, they turn out to be negligible 
in comparison with the other components for all of the situations we have explored and 
therefore in the following we shall focus only on the remaining four larger contributions 
to the current shown in the figures. 

Second, in all calculations we use relativistic kinematics and the full pion propagator 





(58) 
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even when computing the traditional matrix element. The importance of using relativistic 
kinematics is crucial in the evaluation of the angle between h and q arising from the on- 
shell condition given in Eq. fl56|). Non-relativistic kinematics would lead instead to the 
relationship 2h q = Im^io — q 2 , where a factor a; 2 , which is clearly important for the high 
values of q considered here, does not appear. In fact, as pointed in Refs. |], in discussing 
the one-body responses, our approximation to the relativistic current is accurate only if 
the proper relativistic kinematics are used in computing the energy of the particle with 
momentum p. In fact the plane waves are then solutions of the free Klein-Gordon equation 
rather than of the Schrodinger equation, thus automatically accounting for relativity 
in the kinetic energy operator. The additional relativistic dynamics incorporated here, 
arising from the Dirac spinology, enter as modifications of the current operator. For high 
momentum transfers, both ingredients (relativistic kinematics and current corrections due 
to spinology) are of the same level of importance. 

For moderate momentum transfers, say q = 500 MeV/c, the relativistic kinematics 
alone allow one to obtain agreement between the relativistic and traditional K^' TNR 
functions shown in Fig. 2 (compare solid with dot-dashed lines). Although the agreement 
becomes better for approximations NR1 and NR2, all of the curves turn out to be close 
enough to each other to allow the conclusion that for moderate values of q it is sufficient 
to use the traditional seagull current (but including relativistic kinematics) for computing 
one-particle knock-out matrix elements. 

The situation changes at higher g-values. As shown in Figs. 3 and 4, the traditional 
approximation (dot-dashed line) for the time and longitudinal components, KeK^ and 
Reiff, does not agree with the exact relativistic result (solid line). This reflects the 
approximations made in the derivation of this current, where in particular the factors « 
and A are (incorrectly) treated as of higher order. In contrast, our approximations NR1 
and NR2 (dashed and dotted lines) are both very close to the fully-relativistic result 
for all of the g-values considered. Therefore the relativistic corrections included in these 
components of the current in our approximation NR1 (or in its simplified version NR2) 
appear to be sufficient for a proper description of the relativistic effects. 

With regard to the transverse components of the seagull current in the s p = 1/2, 
Sh = — l/2 case, namely KeKf and ImK^ , we first note that for low momentum transfers 
they dominate over the longitudinal components to the left of the quasielastic peak. As 
functions of ui these transverse components are nearly linear and cross the cu-axis some- 
where to the right of the quasielastic peak. This change of sign accounts for the negative 
interference between the seagull and one-body current contributions in the transverse 
electromagnetic response to the right of the quasielastic peak [T5|, |T%fl . The value of u 
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where these functions vanish decreases with increasing q: for q = 2000 MeV/c the posi- 
tion of the zero almost coincides with the center of the quasielastic peak. The figures also 
show that at low energy (below the quasielastic peak) the exact result almost coincides 
with the other three curves. On the other hand, for high energy (above the quasielastic 
peak) discrepancies occur between the exact and traditional currents: thus at the end of 
the allowed u;-region, the traditional current (dot-dashed lines) accounts for only about 
one-half of the exact result (solid lines) in absolute value. On the other hand, our two 
non-relativistic approximations NR1 and NR2 (dashed and dotted lines) are much closer 
to the exact result. Hence our results show that the relativized, simplified approximation 
NR2 to the seagull current is a valid representation of the exact relativistic current for all 
of the values of the momentum transfer considered here. 



3.2. Pion— in— flight current 

Now we perform a similar analysis for the pion-in-flight (or pionic) current j£. First we 
define dimensionless functions K£ for this current as we did for the seagull case; thus in 
the matrix element 

J2(ph'\jr\h'h) = tA^Tx x T 2 ],)F?k* F KZ(q,u,,h) (59) 
h' vm n 

the pionic function K£ reads 

K » ^ W ' h) = 14 fr J Wf { " )M [(P - H') 2 - ml] [(H> - HY - ml] " 

(60) 

We also introduce the traditional non-relativistic function, denoted "P,TNR", for the 
pionic contributions 

(61) 

Finally, our approximated pionic function, introduced in Sect. 2, is simply given by 

K? NR (q, w , h) = -^=K^ R (q, w , h) , (62) 

where "P,NR" stands for non-relativistic approximation for the pionic contributions. 

In Figs. 5-7 we display the various components of the pionic function for the same 
kinematics as employed for the seagull current. The meaning of the curves is the same as 
in Figs. 2-4, except that now only one relativistic approximation is suggested, since the 
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simplicity of our result does not require additional assumptions. The results obtained are 
similar to the ones already found for the seagull current, and we can summarize them in 
the following points: 

• The important contributions arise from the real parts of Kq , K± and K%, and 
imaginary part of K% . The reasons for this dominance are discussed in the Ap- 
pendix, where we are able to inter-relate the general behavior of the curves shown 
in Figs. 5-7 and we explore the symmetries of the various pieces of this current. 

• For low momentum transfers (q below some "moderate" value of about 500 MeV/c) 
the traditional approach, the exact matrix elements and, of course, our present 
approximation are all very close if one takes into account relativistic kinematics. 
This justifies the use of the traditional MEC for low to moderate momentum trans- 
fers ESS S3- 

• For high momentum transfers (q > 1000 MeV/c) the traditional expression (dot- 
dashed line) clearly disagrees with the fully relativistic result (solid line). From 
Figs. 6 and 7 it is apparent that the major differences between these two functions 
disappear if we use our approximated current P,NR (dashed lines). In fact, it is a 
particularly gratifying result that, apart from using relativistic kinematics, the sim- 
ple factor 1 / a/1 + r applied to the traditional current K^' TNR is able to reproduce 
the exact relativistic matrix element remarkably well. Although some disagreements 
between our approximation and the relativistic current for the case of Im exist, 
it is however clear that the traditional result is much worse. Moreover, the dis- 
agreement is found only away from the quasielastic peak where the corresponding 
one-particle emission response is small, the two matrix elements (solid and dashed 
lines) being equal at the peak where the approximations associated with the factor 
1 / \/i + r are expected to work better. 

In conclusion we see that the new pionic current obtained by multiplying the tradi- 
tional non-relativistic one with the spinology factor 1 / a/1 + r significantly improves the 
relativistic content of the current and hence one can use this current for computing one- 
particle emission responses for high momentum transfers within non-relativistic models, 
at least near the quasielastic peak. 

3.3. The large— q limit 

Let us end this section with a brief discussion of the behavior of the currents in the 
large-g limit. We start from the non-relativistic reductions of the currents as given in 
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Eqs. and consider the limit k — > oo. For these conditions 

1 



A 



(2k + r7) 2 + 1 - Jr) 2 + 1 



and 



r = (k + A)(k — A) ~ k(1 — k ■ rj) . 



(63) 



(64) 



By inserting Eqs. (|63"| , |6"4|) into the currents we obtain then for the single-nucleon current 



j" 



1 f G 

G E (l + -k-rj)+i [G M - ) (k x 77) • <x 



— J°k ~ J°k 



J 1 ~ zv^G A f |(cr x k) - ^ [(cr • k)(k x ry) + (cr x k)(k • ry)] 

-|(»7-(«-»7)«)} 

and likewise, for the meson-exchange currents, 
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(71) 



where the inverse of the pion propagator has been expanded to leading order in the 
parameters i] and 
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(72) 
(73) 



At large g, it is of significance that: 



All of the currents grow asymptotically as a/k. This result is supported by our 
numerical results, which show that the currents at q — 2000 MeV/c are roughly 
twice as large as those at q = 500 MeV/c. 



21 



• The time components of the seagull and pion-in-flight currents tend to cancel each 
other. The same happens for the longitudinal components. Hence only the trans- 
verse components of the seagull current survives in this limit. 

• The longitudinal and time components of both the seagull and pion-in-flight currents 
become equal as q — > oo. Since in this limit A ~ k, it follows that in the large-q 
limit the seagull and pionic currents are separately gauge invariant. Moreover, in 
this kinematical regime, the correlations among nucleons are not expected to play 
a significant role, thus implying the separate realization of gauge invariance in each 
sector of the nuclear response. By extension, the current carried by each individual 
meson should be expected to be separately conserved. 

• Finally, if the form factors are neglected, the single-nucleon current and the MEC 
display the same asymptotic behavior in q. Of course the inclusion of form factors 
will change the g-dependences of the currents. 

4. Summary and Conclusions 

In this work we have found new approximations to pionic electromagnetic meson-exchange 
currents using an approach which parallels recent work involving expansions of the elec- 
troweak single-nucleon current in powers of the momentum of the initial bound nucleon 
i] = p/m N . Our goal here and in that previous work has been to obtain current operators 
that can be implemented in computing response functions for high momentum transfers 
in quasielastic kinematics using non-relativistic models. Our approach allows features 
of relativity to be taken into account through the use of relativistic kinematics and the 
Dirac-spinology content implicit in the new currents. 

In this paper we have first illustrated the basic procedure by reviewing the simpler case 
involving the expansion of the single-nucleon current. We have then turned to our main 
focus in the present work and applied the expansion ideas to a study of the pion-exchange 
seagull and pion-in-flight MEC. A distinguishing feature of the recently-obtained single- 
nucleon current is that it incorporates relativistic effects through multiplicative factors 
involving the dimensionless variables k, t and y/1 + r (arising from the Dirac spinology) 
- factors that are easy to implement in traditional non-relativistic models. Accordingly, 
in our expansion of the MEC we have sought to identify corresponding factors which 
can embody the essential features of the relativistic MEC. We have also examined the 
behavior of the currents in the asymptotic limit where q — > oo and found that only the 
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seagull transverse current survives, and moreover that the latter and the pionic current 
are separately conserved. 

Finally, we have tested the quality of our approximations by computing the transition 
matrix elements (ph~ 1 \J^ IEC \F) for the various components of the currents, i.e., for ma- 
trix elements taken between the ground state of a Fermi gas and a particle-hole excitation. 
We have compared the exact relativistic matrix elements with our non-relativistic approx- 
imations and with the traditional non-relativistic expressions. The differences between 
our newly-obtained currents and the exact ones are small even for very high momentum 
transfers, whereas the traditional expressions fail at high q. Due to the quality of our 
results, we believe that these currents can very safely be used in non-relativistic models 
for computing MEC effects in one-particle emission nuclear responses. 

Appendix. Symmetries and relevance of the various 
components of the currents 

In this appendix we study the structure of the integrals K° involved in the calculations of 
the MEC matrix elements to assess the relative weight of the different pieces into which 
they may be decomposed for a variety of kinematical conditions. In fact, in Figs. 2-7 we 
have only shown the dominant among the four components K° (real or imaginary part). 
The components not shown in the figures have been also computed and found to be small; 
here we present arguments to help in understanding why they are so. 

For this purpose it is sufficient to consider only the traditional non-relativistic currents, 
since they are modified simply by multiplicative factors in our various approximations and 
thus do not change our conclusions as far as the issue at stake here is concerned. For 
illustration, we restrict ourselves to the analysis of the pion- in-flight current; the same 
arguments can be applied to an analysis of the seagull current as well. 

Our analysis basically amounts to studying the behavior of the pionic function defined 
by Eq. (fUl). To this end, we first extract the uninteresting constant factor l/k F (2ir) 3 and 
define the following integral appearing in the non-relativistic pionic current 

K» = £ / *h'{P + H- 2HY X Y p ' %^ XSh £^ ' , (74) 

l£J [{P- H') 2 -ml\[{H '- H) 2 -ml] 

which depends upon the dimensionless momenta k, rj and rjp- We aim to identify the 
leading order in the expansion in rj and 1]f of the terms into which this integral splits. 
Since r\ < tjf, we shall identify 0{rj) = 0(t)f)- From this investigation we shall see how 
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the occurrence of cancellations rendering some matrix elements smaller than expected 
emerges. 

We start by performing the summation over the spin of the intermediate hole using 
the completeness relation J2s h , Xs h ,x\ t , — 1 an d by defining the spin- matrix 

r = ^<T.(p-h') Xsh ,xt h ,<T-(ti-h) 

s h> 

= q-(h'-h)-(h-h'f+%x(h'-h)].(T. (75) 

We shall work in the coordinate system where q = qe 3 and study the spin components 
K^i and K^ 2 for which 

T n = q • (h' - h) - (h - h') 2 = q(h' - h) 3 - (h - h') 2 (76) 
r 12 = q(ti -h) x -iq(h' -h) 2 ■ (77) 

Spin component Iffi. In this case we obtain 

K u = I <f*'(P + H - 2HT „ n ^- k }-^,-ll _„ • (78) 



[(P - H') 2 - ml] \(H> - H) 2 - m% 
As a consequence, at the non-relativistic level, we have 

lm#ft = 0. (79) 

This is no longer true for the relativistic pionic contribution, although the corresponding 
imaginary parts of the relativistic pionic matrix elements remain very small. 
Next we examine the real components for /i = 0, . . . , 3. 

KeK^ component 

First, the longitudinal component turns out to read 

q(h'-h) 3 -(h-h') 2 



Re Kl x = J d 3 h'{q + 2h 3 - 2h' 3 



3J [{P - H') 2 - ml) \{H' - H) 2 - m 2 } 
= 0(rj F ) + 0(rj 2 ) + 0{rf F ) , (80) 

where 0(j]p) means that the associated integrand is linear in h or hi . Similarly 0{rf F ) 
means that the integrand is proportional to h 2 , (h') 2 or hh! and so on. The leading terms 
in Eq.(^) are clearly of order rjp and rjp and their expressions are given by 
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0{n F ) = J d 3 ti 



{{P - H') 2 - m 2 w ] [(H> - H) 2 - ml] 
2q{h 3 - h' 3 ) 2 + q{h - h^ 2 
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In comparing these two pieces, we see that the integrand in the term 0(j)f) is proportional 
to h' 3 — hs, which can be positive or negative, thus potentially leading to cancellations. In 
fact, this term is found to be close to zero near the quasielastic peak (QEP). The reason 
is the following: at the QEP p — q, and, for q large, h is almost perpendicular to q and 
thus for the case <fih = 0° at the QEP we can set h ~ he\. The corresponding integral 
accordingly vanishes to first order, namely 

J3W ( h 3 ~ ^3) _ f ,3./ h> 3 



d 6 ti — — / d A ti — = 

[(P _ R' f - ml] [{H> - Hf -ml]~ J [Q 2 - ml] [{H> - H) 2 - ml] 

(83) 

because the denominator is invariant with respect to the inversion of h' 3 . As a consequence, 
for the longitudinal component of the pionic current Ku, we have 

Re Kly ~ near the QEP, (84) 

as we have checked numerically. For instance, in Fig. 6 we see that ReKfy crosses the 
a>axis somewhat short of the middle of the ^-allowed domain (the approximate position 
of the QEP) because of the approximation made in the denominator and of the presence 
of the other term 0(rf F ). This is found to be negative in this region and not entirely 
negligible; hence the total matrix element reaches zero slightly to the left of the QEP. 

Therefore, we see in this case that, although a priori this current if of 0(r] F ), its actual 
weight depends upon the value of the coefficient that multiplies r] F in the expansion, which 
is ^-dependent and in some cases may be small. 

On the other hand, we see in the same figure that in the regions far from the QEP 
(especially for high a;) this component is large: indeed here the cancellations are much 
weaker and then the behavior of ReKfy is again of 0(r] F ). 

Re-K^i component 

The same conclusions obtained for the third component are valid as well for the Re K^. 
In this case: 

Re K° I M L + 2^- - ^ ~ k * ] ~ (h ~ h/)2 

11 ~ J { + 2m N 2m N ) [{P-H'f-ml\\{H'-Hf-ml] 

= 0( VF ) + 0(t? f ) + 0(t? f ) , (85) 



with 



OM ^ /^ [(p _ g , )2 ^^ g)2 _ m;] (86) 



25 



Now, for the same reasons as before, the term 0{t}f) almost vanishes near the QEP. Then 

Re ~ near the QEP. (88) 

Hence Rei^ behaves like ReKf x , as also found in our calculations (see Fig. 6). Even 
finer details of the results can be interpreted in the same manner. For instance: 

• The zero of K® x occurs slightly to the right of the zero of K X1 , because, in the present 
case, a piece proportional to 2q(h 3 — h 3 ) is missing in the term 0(rf F ). Hence the 
0{rf F ) term is less negative than in the K\ x case. 

• Also, to second order, the following relation 

^3 _ <?^o _ [« h , Mh ~ h' 3 ) 2 

k xi --K n jdfi-^pz H , )2 _ m2][{H , _ H)2 _ m2] i»yj 

is seen to hold. For instance, for q = 1000 MeV/c, <ph = 0° and w = 550 MeV, it 
turns out that K X1 = —15, -K X1 = —27, while Kf x = —40. Hence an estimate of 
— 13 follows for the magnitude of the second-order term, represented by the integral 
on the right-hand side of Eq. (p9[) . 



ReK^, ReKf ± components 



Concerning the transverse part we have 

ReK fl = fM2<* - ^ | (p _ff:^:g_ m; | = 0(,f F ) + 0(,f F) . (90) 
The major contribution is expected to arise from the second-order term, namely 

°<*> = / <i3ft ' 2(h - h V-/^-S- ' (91) 

Now, for (ph = 0° and near the QEP, we have h ~ hxex and for the components 1 and 2 
we obtain 

ReK 1 ~ [d 3 ti Hhi-K)h' 3 

11 " J dh [Q*-Tri*)[{H'-H)*-m*}- U {J) 



R«*?i - -/**' ™ —.9.1 ~ ^ " 03) 



[Q 2 - ml\[{H> - Hf - m 2 ] 

In the first case exact cancellations occur when h' 3 — > —h' 3 , but in the second case we 
actually have double cancellations when h 3 — > — h 3 and h' 2 — > —h' 2 . Accordingly it should 
be expected that: 
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• One has 

| Re | <C iRei^nl . (94) 

Actually the cancellations in Reif^ are so strong that this component is even 
smaller than Imif^ for the relativistic current (hence Reif^ is not displayed). 

• Both Re K\y and Re K\ x vanish around the QEP. 

• If <fi h = 90° then h ~ h 2 e 2 : accordingly the roles of ReK 1 and ReK 2 switch, i.e. 

| Re if { x I < | Re | . (95) 

All of these properties have been checked in our numerical results. 

Since Re if^ is zero around the QEP for <p h = 0°, we infer that this piece is actually of 
0(7]%). Moreover, as the approach to zero of Reif? x is faster than in the case of Reif^, 
Re if^ is likely of 0(r)p) (very small). Clearly a precise determination of the actual order 
would require a more detailed analysis of the integrals, or to compute analytically the 
integrals in the static limit. However, we believe that the arguments given above are 
enough for reaching an adequate understanding of the results. 

Spin component iff 2 . For the s p = 1/2, s h = — 1/2 component we have 

iff 2 = [ d*h>(P + H- 2HY- q{h [- h) \;' q{h '- h } 91 . (96) 
12 J V ; [{P-H'f -ml\[{H'-H) 2 -ml] V ; 

Hence the real and imaginary parts of this matrix element are given by 

Re*& = / M(P + H - 2H T [(p _ Rl? - *)._ ^ _ ^ (97) 
Im Kf, = / M(P + H - 2H r w - w ^^ hw _ . (98) 

if° 2 component 

For the time component we have 

K o f ,3 . , / , 2 Jl_ _ 2 J^_\ Q(h> - jQi - iq{h' - jQg 

12 J \ 2m N 2m N ) \(P - H'f - ml}\(H> - Hf - m 2 } 

= 0( VF ) + 0{rf F ) . (99) 



To first order we obtain 

K° ~ />//:: q{h! -h^-iqjh' -h) 2 

K ^-J dh U \(P-H>V-ml\\(H>-HV-ml\ ' (1 ° 0) 
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For <ph = 0° and near the QEP again h ~ h±ei : and hence 

K ° „ [tf h > u q{h> -h) x -iqh> 2 

Al2 " J d kU [Q2 _ ml] [{H> - H) 2 - ml] ' 

As in the case of the K u component, cancellations occur in a way that yields 

Imiffj ~ at the QEP. (102) 
Hence Re K® 2 is expected to be the dominant component. 

K\ 2 component 

A similar result obtains for the longitudinal component, namely 

K\ 2 ~ J d 3 h'q [(p _ q ^ /)2 - _ H) ] 2 _ (103) 

and, for the same reasons as before, 

Im4~0 at the QEP. (104) 
Also the following relationship 

K\ 2 ~ iff* (105) 

is found to hold to leading order, as can be verified in the figures. For instance, from 
Fig. 6, for q = 1000 MeV/c and <ph = 0°, at the maximum (iv ~ 450 MeV) we get 
Re K% 2 ~ -11, ReKf 2 ~ -25 and -K\ 2 = -24.4. 



K\ 2 , K\ 2 components 



For the transverse components we have 



" /^-^ [(P-^iS-^-.i] (1 ° 6) 



and for the real and imaginary parts of K 1 



[(P - H') 2 - m 2 } [{H> - H) 2 - m 2 } 
' : ^ ^'^'[(P - P') 2 - inl}[(H' - H)' 2 - ml] 



7 P-H') 2 -ml } H'-H) 2 -ml} 1 ; 



28 



First we see that the integrand in ~ReK\ 2 is proportional to 2(h[ — h\) 2 > 0, and hence 
no cancellations occur. Accordingly this integral, although of 0(i]p), turns out to be of 
the same order of magnitude as K\ 2 (of 0(i]f), but with cancellations), as is clearly seen 
in our results. On the other hand, cancellations do occur in lmKl 2 ; hence, as before, 

lmKl 2 ~0 at the QEP. (110) 

We thus expect YleK\ 2 to be the dominant part. 
In connection with K\ 2 we obtain 

12 " ~l dh [{P-H'Y-ml\[{H'-HY-ml\ (1U) 

Im^ 2 - fd'ti 2q(h ' 2 - h)2 (112) 

lmA 12 - J dh [{p _ H ,y_ m2][{H ,_ H)2 _ m2] ■ {LU) 

First we see that, as found in our calculations, the following relation 

ReK 2 2 = -\mK{ 2 (113) 

holds, and hence the real part of K\ 2 is negligible with respect to the imaginary part, 
which is thus the dominant one. Also, since the integrand for \m.K\ 2 is proportional to 
ih,2 — h' 2 ) 2 > 0, no direct relationship between lmKf 2 and ~ReK\ 2 exists. 

Again, for (p h = 90°, the roles of K 1 and K 2 are switched, and from the above expres- 
sions we find 

ReK\ 2 {<p h = W) = -Imi^O^gO ) (114) 
Imi^y^O ) = -Rei^^gO ) . (115) 
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Figure 1: MEC diagrams considered in this work. 



500 MeV/c, (f> h = 0° 
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Figure 2: Seagull current matrix element K s ^ for q = 500 MeV/c. The kin- 
ematics for the hole h is h = 175 MeV j c and <f>h = 0°. First column: spin 
{\,\) component. Second column: spin {\,—\) component. Solid: relatwistic, 
dashed: our approximation NR1, dot-dashed: traditional non-relatwistic TNR, 
dotted: our approximation NR2. 




Figure 3: The same as Fig. 2 for q 



= 1000 MeV/c 



q = 500 MeV/c, 4> h = 0' 
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Figure 5: Pion-m-flight current matrix element K v for q = 500 MeV/c. The 
kinematics for the hole h is h = 175 MeV j c and <f> h = 0°. First column: spin 
(\,\) component. Second column: spin {\,—\) component. Solid: relatwistic, 
dashed: our approximation NR, dot-dashed: traditional non-relatwistic TNR. 



